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We propose a new application of the Gutzwiller trace formula formalism, to give a compact 
\ expression for the semiclassical vacuum pair production rate in quantum electrodynamics, 

, for general inhomogeneous electromagnetic background fields. 
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■ The Gutzwiller trace formula has found a wide variety of applications in theoretical and math- 

ematical physics . Here we point out a new area where its language provides insight and 

simplification to a difficult computational problem in relativistic quantum field theory. Vacuum 
D \ polarization effects in quantum electrodynamics (QED) predict that electron-positron pairs can 



be produced from the vacuum in the presence of a classical electric fie^ 
nomenon was predicted, and its rate estimated, for uniform fields in [4, 



. This remarkable phe- 



u. im 

,3,3, 



but has not been 



directly observed as the rate is tiny for accessible field strengths. It is conceivable that sufficiently 



strong electric fields may be reached in X-Ray Free Electron Lasers 



71], but such fields have strong 



> 

o 
o 

' temporal and spatial inhomogeneities. Unfortunately, very little is known about the rate when 

o : 

. the background field has such general inhomogeneities. On the other hand, in the approximation 

O 

where the background electric field has a fixed direction and a magnitude that varies in just one 
^ . dimension, either spatial or temporal, one can use WKB-based techniques 8,3, [l^. A promising 

\ approach for going beyond this one-dimensional case is the "worldline instanton" method [ill, [l^. 

)ased on an instanton approximation to Feynman's worldline path integral formulation of QED 
1^ . Another related approach is a direct Monte Carlo evaluation of the worldline form of the 
effective action |14|]. In this note we propose a new approach to this problem, based on a close 
connection between the worldline instanton approach and the Gutzwiller trace formula 
This connection gives a well-defined computational strategy for treating multi-dimensional inho- 
mogeneities in the background electromagnetic field. 

The technical problem is to compute the imaginary part of the effective action in the clas- 
sical electromagnetic background field, from which the vacuum pair production rate follows Q]: 

1 _-2Tmr nT._„-n 1_ r £ _1 J _r c . 

production 



-^production = 1 — 6 ^imT 2Imr. For example, for a constant electric field of magnitude £, the 



leading weak field result (we consider scalar QED) is 



Vol 
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The basis of our proposal is the worldline formalism of QED |l3l . [isl . [igI , in which the effective action 
is expressed in terms of a quantum mechanical path integral in four-dimensional Euclidean space, 
with paths x^{t) parametrized by proper-time r. This approach has led to many beautiful advances 
in our understanding of perturbative scattering amplitudes [3], but here we propose to use it to 
extract non-perturbative information. The effective action for a scalar charged particle (charge e, 
mass m) in a Euclidean classical gauge background A^{x) is the functional {D^ = + ieA^ is the 
covariant derivative): 



r[^] 



tr In {-D^ + rr?) 



(2) 



T 
dT 



e-n-Dl) |^(0)^ 



""'^ I d^x^o) J Vx exp 



T 



dT\^ + e A^x^ 



x{T)=x{0)=x(0) 

In the last line, the trace of the associated Euclidean propagation operator has been written 
as a functional integral J Vx over all closed Euclidean spacetime paths x^(t) that are periodic 
(with period T) in the proper-time parameter r [3]. We use the QED worldline path integral 
normalization conventions of [191 . 

The strategy of the worldline instanton method [U] is to evaluate the quantum mechanical path 
integral in ([2]) semiclassically and then to evaluate each of the T and x^'^^ integrals by steepest 
descents. These are precisely the steps in deriving the Gutzwiller trace formula , although 

there one is concerned with a non-relativistic Schrodinger operator rather than the Euclidean Klein- 
Gordon operator, an oscillatory amplitude e^^^^ rather than the Euclidean form e"*^, and the trace 
of the resolvent rather than the trace of the logarithm. Nevertheless, despite these differences, in 
this note we show that the worldline instanton computation can usefully be formulated in the 
language of the Gutzwiller trace formula. 

The first step is to make a semiclassical approximation for the propagation kernel 



K{x,x';T) := (x|e-^(-^M) 



det 



_cPB_ 
dx dx' 



-R{x,x';T) 



(3) 



where R{x, x'; T) is the Hamilton principal function for the classical trajectory from x to x' in four- 
dimensional Euclidean space, in the proper-time interval T. This classical trajectory is obtained 



3 



by solving the Euclidean classical equations of motion 

Xf, = 2eF^^{x) , {fi,u=1...4) , (4) 

where F^i, = d^A^ — dyAfj_ is the background field strength. To evaluate the trace in ([2]) we will need 
the diagonal propagation kernel K{x^'^\x^'^'> ]T), but for now we consider the point-split propagation 
from X to x' . The classical equations of motion (jl]) are those for a charged particle moving in an 
inhomogeneous electromagnetic field F^^{x), so the "energy" is conserved on a classical trajectory 
: E = jxjj^ = constant. 

The next step is to perform the T integral by steepest descents. The critical point of the 



exponential factor arises when ^ 



-m 



This has a natural classical interpretation in terms 



of the Legendre transformation between the Hamilton principal function R{x,x';T) [expressed in 
terms of the total time elapsed along the trajectory] and the action S{x, x'; E) [expressed in terms 
of the constant energy of the trajectory] : R{x^ x';T) = S{x, x';E) — ET. It follows that |y = —E, 
and ^ = T. Thus, the critical point Tc of the T integral occurs when E = vn?, so that 

1 



l^^^e--'^K{x,x';T) 



(2^)2 



det 



dx dx' 



27r 
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pie way if we consider coordinates x'f'^ along the classical trajectory, and x^'f^ transverse to the 



up to a possible phase that we discuss later. The two prefactor contributions combine in a sim- 
nsider a 

trajectory. Then [l|, 3, Q] 

\ dxj_ dx\ 



det 



/ d'^S{x, x'] m?] 



Xll X' 



(6) 



The final step is the coincident limit x — > x' = x^^\ and trace over x^^\ This trace is also 
done by steepest descents and implies that the closed loop is in fact periodic . Periodic 

solutions to @ are known as worldline instantons [l2]. From ([6]), the integration over Xij*^^ yields 
a factor J = Tc/2 [reparametrization invariance of the periodic orbit], while the x^^ 
integral produces another determinant factor. Remarkably, this determinant factor combines with 
the remaining transversal determinant factor in ^ to give: 



det 



( d'^S{x,x';m?) _|_ d'^S[x,x';nr') _^ d'^ S{x,x' ■,m?) _|_ d'^ S[x,x' -^m?) 



\ dx^ dxj^ 



dx'. dx± 



+ 



dx I dx', 



+ 



dx', dx', 



det f 

\^ dxj^ dx'j^ 



det 




--: det (1 - J) 



(7) 



where all determinants are evaluated at vanishing transverse displacements. Here J is the mon- 
odromy matrix, for a 6-dimensional surface of section in phase space transverse to the periodic phase 
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9 / '^^i 

space orbit with constant energy E = m . Consider an initial transverse displacement I 

w± 

from a point on the closed orbit in phase space, and evolve for time T, and the final displace- 

/ 6x'\\ ( 8x'^ 

ment from the orbit is related to the initial one by the monodromy matrix: I I ~ I 

Putting all these parts together, and collecting phases carefully [ij], one obtains a compact final 
expression : 

IniTp^ ^= . (8) 

V^det (1 - J) 

The principal advantage of expressing the computation in this language of the Gutzwiller trace 
formula is that the total prefactor is encapsulated in a single determinant, which moreover has a 
natural mathematical and geometrical meaning in the Euclidean phase space. In previous work 
9. [lol. 1^ the various prefactor contributions have been evaluated separately, and then combined 
at the end. Thus, the computational strategy is as follows: 

1. solve the classical equations of motion in four dimensional Euclidean space to find all closed 
periodic trajectories of energy E = vn? : the "worldline instanton(s)" . 

2. evaluate the classical action S{E = m?) on these trajectories. The dominant contribution 
comes from the trajectory(ies) with largest _ 

3. compute the prefactor from the monodromy matrix J for the dominant trajectory (ies). 

m r for the case of a one-dimensional 

inhomogeneity, which can be computed in several other ways 
case of a time dependent electric field directed in the X3 direction. We can choose a Euclidean gauge 
field A^{xii) = § /("^ ^4)) where £ characterizes the overall magnitude of the associated electric 
field, 00 characterizes the scale of the time dependence, and /(a; 2:4) is some smooth function. For 
example, for a constant electric field £{t) = we have f{x) = x; for a sinusoidal electric field 
£{t) = £ cos(a; t), we have f{x) = sinh(x); and for a single-pulse electric field £{t) = iSsech^(a; t), 
we have /(x) = tan(x). Then the classical action on a periodic trajectory of energy E can be 
written [here y := /(x)] 



The only concrete comparison we can make is to compute 



(e£ „, _ /-^ 



,2 



This is precisely the exponent appearing in the standard result for the pair production rate 
Q,Q. To evaluate the prefactor, we can choose Then the transverse xz direction is in 
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fact an invariant "flat" direction, so we do not need to perform the transverse integration. This 
illustrates the important point that ([8]) must be interpreted appropriately when there are physical 
zero modes. Thus, we go back to ([5]) and observe that = — Furthermore, the other 
determinant factor in ([5]) is easily computed (see [l3]b) using the Gel'fand-Yaglom formula: 



det 



dx dx' 



m 



1 



16 r2 



Xa 



4 \dE^ 



(10) 



Thus, relative to the constant spatial volume V3, 



ImF 



V3 2(47r)2m 



-S{E) 



95 d'^S 
dE\J dE^ 



(11) 



Note that (jlip agrees precisely with the conventional WKB result 

We now turn to a multi-dimensional example. Consider the two-dimensional Euclidean problem 



Q in the X3-X4 plane, with F34 = F{r), where r := \/x'l + x\. The associated Minkowski electric 



field points along the axis and is a function of y^t^ _ x'i^, i.e. a configuration studied, e.g., in [isl ]. 
There can exist circular orbits centered around r = 0. The fiuctuation determinant for fiuctuations 
around such an orbit of radius rg follows from the corresponding monodromy matrix J. In polar 
coordinates the circular orbit is characterized by r(T) = tq, and 9[t) = 9q = 2Fq := 2F(ro). 
Linearizing the equations of motion in fiuctuations p and around the periodic trajectory, where 
r(r) =: tq -|- p{t) and 6{t) =: Oq{t) + '(^(t), and solving the resulting equations for the intial 
conditions p(0) = 6x'j_ and p(0) = 5p'j_, leads to the following solution for the radial fiuctuations: 



fe'_L cos(2rFoo-) + 6p_^{2Foa)-^ sin(2rFoa) 



(12) 



Here o" := [1 -|- ^{dr\roF)]^/'^ , and we made use of p9o = 'Ovq, which follows from the conservation 
of the magnitude of the velocity \/ x^ + x\. To compute the transverse deviation from the orbit 
after one cycle, in principle, we have to calculate the time needed in order to return to the same 
longitudinal coordinate {5x\\ = 0), which here means 6{T + 5T) := 27r, where T = tt/Fq is the 
period of the unperturbed orbit. Putting into Eq. (jl2p r = T + 6T, instead of r = T, however, 
leads merely to corrections quadratic in the initial fiuctuation parameters 5x'j_ and 6p'j_. Therefore, 
6q'^ = p{T), and 6p'[^ = p{T). Then the monodromy matrix is: 



J 



cos(27r(T) 2^00" sin(27r(T) 

-(2Fo(T)-isin(27ro-) cos(27ro-) 



(13) 



The corresponding fiuctuation determinant is given by det(l — J) = 4sin2(7r(7). 
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We conclude with some comments and open problems, (i) Finding closed periodic orbits to (jH) 
is non-trivial, but recasting the problem in phase space proves helpful, (ii) If the physical electric 
field is too localized in space, then we know physically that the pair production rate vanishes 
[since the virtual vacuum dipole pairs cannot gain enough energy from the field to become real 
electron-positron pairs]. In simple cases this corresponds to the non-existence of periodic classical 
Euclidean trajectories [3]. It would be interesting if this were more generally true: that the mere 
existence of such worldline instanton loops might be used as an indicator of pair production, (iii) 
The phases arising from the steepest descent integrals combine to give Im F in the cases where the 
electric field is a function of either t or x {l^ . but the general mixed case needs further analysis, 
(iv) If the gauge field corresponding to the external field can be put into the non-linear gauge 
where A'j^{x) = constant (= E), then we can solve the simpler first-order equations = —2eA^{x), 
as was observed long ago by Nambu [l^ . (v) It would be interesting to extend our method to 
inhomogeneous non-abelian fields, for which little is known beyond simple quasi-abelian cases. 
This suggests studying the Wong equations [2l describing the classical motion of a color-charged 
particle in a non-abelian background, which is a much richer mathematical system. 
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